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Abstract 

We calculate masses, production cross sections, and decay rates in the 
Two Higgs Doublet Model of type II. We also discuss running coupling 
constants and Grand Unification. The most interesting production 
. channels are gg — > hP,H°,A° on mass shell, and qq,gg — > h°Z and 

O ! QQ' ~ > hPW^ in the continuum (tho there may be peaks at m^o). The 

most interesting decays are h°,H°,A° — > 66-jets and t + t~, and, if 
above threshold, H° -> ZZ, W + W~ and h°h°. The following final 
states should be compared with the Standard Model cross section: 
bbZ, bbW ± , t+t-Z, t + t-W ± , bb, t + t~, ZZ, W + W~, 3 and 4 b- 



jets, 2r + +2r - , bbr+r', ZW + W~ , 3Z, ZZW ± and_3W ± . Mass 
r~| , peaks should be searched in the following channels: Zbb, ZZ, ZZZ, 

bb, 46-jets and, just in case, Z7. 
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Among the extensions of the Standard Model that respect its principles and 
symmetries, and are compatible with present data within a region of param- 
eter space, and are of interest at the large particle colliders, is the addition 
of a second doublet of higgs fields. In this article we consider the Two Higgs 
Doublet Model of type IIpQ. The higgs sector of the Minimal Supersymmet- 
ric Standard Model (MSSM) is of this type (tho the model of type II does not 
require Supersymmetry) . The physical spectrum of the model contains five 
higgs bosons: one pseudoscalar A° (CP-odd scalar), two neutral scalars H° 
and h° (CP-even scalars), and two charged scalars H + and H~ . The masses 
of the charged Higgs bosons m H , and the ratio of the vacuum expectation 
values of the two neutral components of the Higgs doublets, tan (3 > 0, are 
free parameters of the theory. 

In [2] we obtained limits in the (ra^, tan/3) plane using measured de- 
cay rates, mixing and CP violation of mesons. In this article we present 
graphically the corresponding limits on m#o, m^o and m^o. Then we cal- 
culate production cross sections, decay rates and branching fractions of the 
higgs particles. Next, we obtain the running coupling constants and discuss 
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Grand Unification. Finally, in the Conclusions, we list interesting discovery 
channels. 



2 Masses 

The masses of the neutral higgs particles as a function of the masses of the 
charged higgs m#, tan f3 and the masses of Z and W, calculated at tree level, 
are: 

2 2 2 /i \ 

m A0 =m H - m Wl (1) 




o 2 2 2,2 

2m h0 = m H -m w + mr z 



(2 2, 2\ 2 a 2 / 2 2 \ /tan 2 /5 l\ 



(3) 



We have re derived these equations in agreement with the literature. [1 

In reference j2] we obtained the limits in the (tan (3, ran) plane shown in 
Figure H From that figure and Equations 0121 and El we obtain the limits on 
the masses of the neutral higgs particles shown in Figure 121 

Radiative corrections can be very large. In the MSSM the largest contri- 
butions arise from the incomplete cancellation between top and stop loops. 
The corresponding plot similar to Figure El with radiative corrections can be 
found in 0. 

3 Feynman rules 

The Lagrangian for the VHH interaction is: pQ 
1 9 



Cvhh = ~Y-W+-H~ d» [H° sm{a - (3) + h° cos{a - (3) + iA°] + h.c. 



. -Z„{iA° d " [H° sin(a - (3) + h° cos(a - (3)] 
2 cos uw 

- (2 sin 2 6 w -l) ■ H- <r d^H + } (4) 
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Figure 1: Lower and upper limits on tan/3 as a function of the mass of the 
charged higgs m# from meson decay, mixing and CP violation (continuous 
curve) compared to limits obtained by CDF [3 , D0[4J and LEP2[5 , all at 
95% confidence. Taken from 
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Figure 2: Allowed regions of the masses of the neutral higgs h°, H° and A 
as a function of the mass ran of the charged higgs H^. From Figure Hand 
the tree level Equations [TJ El and El Radiative corrections raise the allowed 
region of h°.6\ 
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where 

Ad^B = A(d^B) - (d> x A)B. (5) 
The Lagrangian for the VVH interaction is: 

Cvvh = (gm w W+W-» + -^-Z,zA 
\ * 2 cos Ow J 

x [H° cos(/3 -a) + h° sin(/3 - a)] . (6) 

There are no vertices ZH°H°, Zh°h°, ZA°A°, ZW+H' or ZH°h°. The 
interactions of neutral higgs bosons with up and down quarks are given by: 

£>AHhff — 7. — n [ufVf(H° sin a + h° cos a) — icos/3 ■ Uf^vjA ] 
2m\v sin p 

gm r 



2m w cos (3 

x [upVf,(H cos a - h sina) - isin/3 • tifrfvpA ] (7) 

where f — u,c, t, v e , v^, v T and /' = d, s, b, e~ , r~ . The Lagrangian corre- 
sponding to the H ± ff vertex is: 

W = ^t^^ H+v ff'^( A + B ^ v f' 

+H-V; f ,u f ,(A-B 1 5 )v f } (8) 

where A = (my tan (3 + m,f cot ($) and B = (my tan (3 — rrif cot (3). Vff is an 
element of the CKM matrix. The Lagrangian corresponding to three higgs 
bosons is: 



C 3h = —gH{H + H 



mw cos(/9 — a) — — cos(2/9) cos([3 + a) 

2 COS 0\\r 



+ mz n H°H° cos(2a) cos(/3 + a) 
Acos6 w 

+ mz n h°h° [2 sin(2a) sin(/3 + a) - cos(/3 + a) cos(2a)] 
4cost^vK 

mz A°A°cos(2/3)cos(/3 + a)} 



4 cos 6^ 
-gh°{H + H 



m w sin(/9 — a) + ^- — cos(2/9) sin(/3 + a) 

2 COS uyy 



+ mz n h°h° cos(2a) sin(/3 + a) 
4 cos 6>vf 



m z 



4 cos 



[2 sin(2a) cos(/5 + a) + sin(/3 + a) cos(2a)] 



+-^—A°A°cos(2(3) sin(/5 + a)}. (9) 
4 cos &w 




Figure 3: Feynman diagrams corresponding to the production of h° in the 
channel qq — > h°Z°. 

Vertexes with four partons including two Higgs bosons are 

U = e*A^H+H- + e9C0S ^ w) A^H+H- 

cos v w 

sin (p - a) A^W^H°H T + ^ cos (/3 - a) AJ¥ ±lx h H* 
±^A lx W± ti A°H*, (10) 
The H + H~ / j vertex is 

Ch+h-^ = ~ig sin e w A tl H- < d^H + (11) 

The higgs propagators are: ij (k 2 — m? + is). 

Feynman diagrams corresponding to the production of Zh° are shown in 
Figures E] and Note that the invariant mass of Zh° can have a resonance 
at m^o which is an interesting experimental signature. Feynman diagrams 
corresponding to the production of W ± h° or W ± H° are shown in Figure |3 
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Figure 4: Feynman diagrams corresponding to the production of h in the 
channel gg — > h°Z°. Continued in Figure El 



4 Decay rates of /i° 

Calculating the Feynman diagrams of Figure U\ we obtain the decay rate 
corresponding to h° — ► gg: 



m^o cos 2 



' a 



647r 3 sin 2 /3 
+r r tan f3 tana [(r T — 1) /(r T ) + 2] 
-r t [(r t -l)/(r t ) + 2]| 2 



|r& tan /Ataxia [(t& — 1) /(t&) + 2] 



where 



4m, 



2 ' 
h° 



in 



(note that n will change from Section to Section), 

I —2 arcsin ( t~ l/z 1 



ln U(i-n)^ 



27T 



if Ti > 1 

if Ti < 1, 



tana 



1+F Y 
l-Ff 



(12) 
(13) 



(14) 



(15) 
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Figure 6: Feynman diagrams corresponding to the production of h°W ± 
H°W ± . 
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VTLr 



m 



w 



rn 



H 



rn 



H 



G 



1 + 



nir 



rn 



w 



rrir 



1 - 



m 



tan 2 (3-1 



m H 171 h J m H V m H J V^an (3 + 1 
Calculating the Feynman diagrams of Figure |H1 we obtain: 



and 



r(h° 



T(h c 



T(h° 



^ ?>G FTn 2 c m h o cos 2 a ^ I/;.', 



2X1 



V2 ■ Air sin 2 (3 

3 

3GFmlm h o sin 2 a / 4m 2 \ 2 
' y/2 ■ 4vr cos 2 /3 V ™J<J 

_^ + , _ G F m 2 m fe o sin 2 a / _ 4m 2 \ 
T ' ~ V2-4n cos 2 p \ m 2 h0 J 
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f 




f 

Figure 8: Feynman diagram corresponding to the decays h° — > bb, cc, rf. 
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5 Branching fractions of h° 

From the preceding decay rates we obtain the following branching fractions 
for the case rrib,m c , m T <C m h o < 120GeV/c 2 : 

B (h° -> 66) = 3m fe sin2 a (21) 

3m 2 sin 2 a + 3m 2 cos 2 a cot 2 (3 + m 2 sin 2 a + J 



and 



m 2 sin 2 a 



B (h° -> r+r~) = — = " (22) 

3m^ sin a + 3m 2 cos 2 a cot /? + m 2 sin a + J 



where 



a 2 m? cos 2 a. „ v-^ \ ( 1 f, /r,\ 1 2 

J = HHr 5-: tan /3 tana > r 7 - — InM +w +2 

8vr 2 tan 2 5 1 H ^ M 2l U/ / 

j=b,r 

-r t {(r t -l)f(r t ) + 2}\ 2 . (23) 

For 90 < m H < 1000GeV/c 2 , B (h° -> 66) varies from 0.856 to 0.944. Ne- 
glecting f? — > gg) and the contribution of cc we obtain 

fi ( h o 6 ft) = ^m 2 = 0.944 (24) 
v 7 3m£ + m 2 



and 



2 

mi 



B (h° -> r+r") = = 0.056. (25) 

' 3m^ + m 2 



6 Decay rates of H ± 

The tree level Feynman diagram of Figure El gives the following decay rate: 

r (If* W** ) = ^ C ° S2a 2 m [l + tan/3tana] 2 
v ; 167rm^(l + tan 2 /5) 1 J 

xA 3 ' 2 (m%,rr&,m 2 h0 ) (26) 

where 

A(o, 6, c)=a 2 + b 2 + c 2 - 2ab - 2bc - 2ca. (27) 
Similarly from the Feynman diagrams of Figures E3 and ^2 we obtain 

r(if ± ^W ± if°) = V2G F (tan P -tana) 2 



lQ7cm 3 H (1 + tan 2 (3) (1 + tan 2 a) 
xA 3/2 (m 2 H ,m 2 w ,rn 2 HQ ) , (21 
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w 




Figure 9: Feynman diagram corresponding to the decay H — > W ± h°. 



W 




Figure 10: Feynman diagram corresponding to the decay H^ 1 —>■ W ± H° . 



r (H± - W±A°) = ^P^A 3/2 Kc <, ml) . (29) 

7 Decays of H°. 

The tree level Feynman diagrams of Figure H2l give the following decay rates: 

8tt V m H oJ 

where Nj = 3 for quarks, Nf = 1 for leptons, = sin 2 a:(l + cot 2 /3) for 
f — u,c, t, and B 2 ? = cos 2 a(l + tan 2 f3) for f — d,s, b, e~ , r". 
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Figure 11: Feynman diagram corresponding to the decay H ± — > W ± A°. 



From the Feynman diagrams of Figure EH we obtain 



y/2G F a 



2 m? 



where 



r(H «^ gg) = """^f "° 3 m 2 a(l + cot 2 ,3) 
|^5> [(T( _ 1)/(T()+21 

i=b,T 

+n[(r t -l)/(r t ) + 2]| 2 (31) 



4m 2 , , 

n = -V". 32 



From the Feynman diagram of Figure E] we obtain 

,2 \ 1/2 n + _2„v2 



a J 

+ tan 2 0) 



327rm H o V / (l + tan 2 a) 3 (l 

^ — (tana + tan/5) — (1 — tan a tan /3) .(33) 

1 — tan a 

From the Feynman diagrams of Figures and we obtain 

32nm H o \ m HQ ) 



°H° 

( tan2 ^~ 1 ) ,[1 -tan a tan/?] 2 , (34) 

^l + tan 2 a)(l + tan 2 /3) 3 

\P1Gf (1 + tan/Jtana) 2 ??^^ 
32tt(1 + tan 2 f3) (1 + tan 2 a) 

x (12a; 2 - Ax + l) (1 - 4a;) 1/2 (35) 



r (h° -> zz) 
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where 



rn 



x = (36) 



Similarly, from the diagram of Figure El we obtain 

t/2Gf (1 + tan /3 tan a) 2 m 3 ^ 



16vr (1 + tan 2 3) (1 + tan 2 a 



x (!2y 2 -% + l)(l-%) 1/2 (37) 



where 



m 



w 



m 2 H0 



From the Feynman diagram of Figure we obtain 

y/2G F (t an a — tan 3) 



From the Feynman diagram of Figure IT9"1 we obtain 



Aitm H o (1 + tan 2 /3) (1 + tan 2 a) 



x 



_ . (1 — tan 3 tana) 1 — tan 2 ' 
1 + tan a tan a) — 



2 cos 2 0^ 1 + tan 2 8 



Am 2 ^ 1/2 



(3? 



167rm^ (1 + tan 2 a) (1 + tan 2 3) 
xA 3 / 2 (m 2 H0 ,m 2 w ,m 2 H ). (39) 



The Feynman diagrams corresponding to h°, H° — > 77 are shown in Fig- 
ure HB 

r(F° -> = 0. (41) 

8 Decay rates of A 

From the tree level Feynman diagram of Figure |^ we obtain 



Vggr— 2 
16vrm3 (1 + tan 2 /3) 

(42) 



r (A - = - V ^ cos a [1 + tan/? tan «] 2 A 3 / 2 (m 2 ,m 2 o,m| 
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f 



Figure 12: Feynman diagram corresponding to the decay H [ 




Figure 13: Feynman diagrams corresponding to the decay H 1 




Figure 14: Feynman diagram corresponding to the decay H { 
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A 



Figure 15: Feynman diagram corresponding to the decay H° — > A A . 




Figure 16: Feynman diagram corresponding to the decay H° — > ZZ. 




Figure 17: Feynman diagram corresponding to the decay H° — > W + W . 



17 




Figure 18: Feynman diagram corresponding to the decay H [ 



H + 



H 



H" 



Figure 19: Feynman diagram corresponding to the decay H l 
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Figure 20: Feynman diagrams corresponding to h°, H* 
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2\ 1/2 
7 



N f 



(43) 



From the tree level diagram of Figure E2 we obtain 

T (A" -> //) = — - — m A om f A f 1 r 

87T \ m lo/ 

where Nf = 3 for quarks, iVy = 1 for leptons, Af — cot /? for / = it, c, £, and 
Af — tan /3 for / = d, s, b, e~, r~ . 

From the Feynman diagrams shown in Figure ESI we obtain: 

V2G F a 2 - 3 / - 2 ^ 3 



T{A°^Z 1 ) = 



em m A° 



5127T 3 sin 2 9w cos 2 6w 

-sin 2 # 
2 3 

(l- 

2 3 



1 - 



rrir 
m 2 A o 



| tan/3 



/l 4 
2 cot /3 - - - sin 2 W 



^ v vv \ A u / 

) / (n, A 6 ) + Q - 2 sin 2 W J 7 (r T , A T ) 



where 



7i = 



4m 2 



A, = 



I(T t ,A t J 

Amf 



m 2 z 



and 



^(ri,A 4 ) = -^-{/(r i )-/(A 4 )} 



A, - Tt 

From the Feynman diagrams of Figure EH we obtain the decay rate 
T(A° -> ^) = ^J^l tan/3 ]T 75/(7*) + cot/3r, f(n)| 2 . 

r (a -> W) = r (a -> #°#°) = o. 

From the Feynman diagram ESI we obtain 

r (a° - w^) = a 3 / 2 Ko,^,0 . 

From the Feynman diagram EH1 we obtain 

T (A -> = v^Gf (tan/3-tana) 2 

^ ~* ' 16vrm3 o (1 + tan 2 a) (1 + tan 2 (5) 
xA 3/2 (m A0 ,m 2 H0 ,m 2 z ) . 
From the Feynman diagrams of EH we obtain: 



(44) 
(45) 
(46) 

(47) 
(48) 

(49) 



(50) 



r(A°^77) = 



V2G F a 



em 



256vr J 
tan /? 



^nf(n) + r T /(r T ) 



cot /? 



(51) 
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Figure 21: Feynman diagram corresponding to the decay A 
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Figure 22: Feynman diagram corresponding to the decay A 




Figure 23: Feynman diagrams corresponding to the decay A 
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Figure 24: Feynman diagrams of A — > gg. 



W 




Figure 25: Feynman diagram of A — > W ± H Zf . 



Z° 




Figure 26: Feynman diagram of A — > 
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Y 

Figure 27: Feynman diagrams for A — > 77. 



9 Decay Z —> h 0r y 



From the Feynman diagrams of Figure 1281 we obtain: 

r (Z -> y/2G F al m m% ( 

K 1] 64vr 3 sin 2 6 W cos 2 6 W \ 



cos 2 a cos 2 j3 



where 



Ti = 



cos (3 sin /? 
+ tana 



647T 3 sin 2 # w cos2 6w z , 

j~ . ^ [tan a tan /? ( - - | sin 2 ^ J F(r 6 , A, 
? sin p \ 2 3 / 

tan /3 Q - 2 sin 2 W J F(t t , A t ) 

-2 Q-lsin 2 ^ F(r t ,A t )] 

r 1 — tan 2 /3 tan (3 + tan a 

tanp — tana H =— — 

H l + tan 2 /3 2cos 2 ^ 

^(l-2sin 2 M^,A H ) 



+ 



— n 

— — (tan (3 — tana) cos 2 9w[4 (3 — tan 2 9w) K( T w, ^w) 
+ { f 1 + ^) tan ' ^ ~ ( 5 + ^) } J(ttv, A w )]| 2 (52) 



Amf 

2 ' 
m 2 



A; 



4m 2 



Ti? = 



4m 2 



2 ' 



Ah = 



4m 2 H 

2 ' 



(53) 
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1 Tiki 



2r i -A i (n-Ai) 2 
Ti - A,- 



+- 



1 TiAi 



4 Ti - A,: 



1 + 



I(t h ,A h ) 



(x — 1) arcsin (2 1 ^ 2 ) 
i(l-a;) 1/2 [in ' 



fo(T,)-0(A,)} 

{/fa) -/(A,)}, 

if x > 1 



l-(l-x) 1 /2 



27T 



if x < 1, 



1 r g A g _ 

2 (th — A#) (r^-A^) 2 



4 , [/(t*) - /(Ah)] 



4 r H - A H 



TW 



Am 2 w 



m 



2 ' 
h° 



A 



ir 



4m^ 



m 



2 ' 



r w A 



w 



[f(r w )-f(A w )}. 



(54) 



(55) 



(56) 

(57) 
(58) 



A(t w - A w ) 

The decay width of Equation EH turns out to be negligible compared to the 
full width of Z so we can not use it to constrain the mass of h°. 



10 Vertex with four particles 

The decay rate corresponding to the Feynman diagram |2H is: 



3(7^ sin 2 9w (1 + tan (3 tan a) 2 m 5 



32 (1 + tan 2 /3) (1 + tan 2 a) tt 3 (a^) 1/2 
x{iAl (l, + < + ^ 



Ah{l, x w,x h °)+l-x™-x h * 



2 («) 



In 



H 


x H 


h° 

x H 


x H 



X 



'if ^H 



2 «^°) 



} (59) 
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Figure 29: Diagram for H ± — > W /± 7/i°. 



where 



11 Production of #° and A 

From the Feynman diagrams in Figure EDI we obtain 

vr 2 r (A -> 2p) 7 A Z" 1 ^ 



CT -> AX) 

4tt 2 7a 



to. 



/ 9 {x a ,m ) g I — ,m J 



+ 



3m A 



,q=u,d,s,c,b 



qq) 



1 eke 



7,1 



•r, 



(x a , m 2 ) /, 



7a 2 
— , m 

•En. 



(61) 



where A = h°,H°,A° and 7a = m 2 A /s. Here / g is the unpolarized parton 
distribution function for quark or anti-quark q and g is the parton distribution 
function for gluons. m 2 is the factorization scale. T(h° — > gg) is given by ()12j) . 

r(#° -> by ||3D, r(A° -> by (HZ|), r -> cc) by (HHD, r -> bl) 

by (O, r(#° -> qq) by (J3UJ), and finally, T (A -> gg) is given by 



12 Production of /i u Z u X 



A production channel with interesting experimental signature is pp — > h°Z°X. 
The differential cross section obtained from the Feynman diagrams in Figure 
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Figure 30: Feynman diagrams for pp — > AX with A = if , A . 
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His 



d 2 a 



dyd (pj 



where / is q or g, 



x§(//^°) 



m^ - k 



x 6 



y/smTe y + m 2 h0 — m 

r~ —'a 
s — y'srriTe y 

m T = (m| + p 2 T ) 2 , 



x a y/srriTe y + m 2 — m 
x a s — -J~smT£ y 



S X a XbS, 



(62) 



(63) 

(64) 

(65) 
(66) 



and 



2 A (s, m 2 , sin 2 9 
Pt= 4l ' 



m 2 ,m|)) 



(67) 
(68) 
(69) 

y is the rapidity, 9 is the angle of dispersion, and pt is the transverse mo- 
mentum of Z°. For the light quarks u, d and s we obtain 



u — - {rn 2 h0 + m| — s — cos6A 1 / 2 (s, 



»; 2 2,-2 

ut = m h0 m z + sp T . 



da 
di 



-Glm\ 



sin 2 (P — a) 



48tts~ f z (s _ ml) 2 + m |r| 

2 



3m| + 



A (s,m 2 ,m|) 



sin^ 



(70) 



where gr^ = t 3 £ (/) and #y = (/) — 2g^ sin 2 9w Coefficients in Equation 
()70)1 are given in Table H The Standard Model cross section is obtained by 
omitting the factor sin 2 (/3 — a) in Equation (|7()J) . The contributions to the 
cross section from the heavy quarks c and b are negligible. is the total 
decay width of the Z°. 
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/ 


h L (/) 


9/ 






e", H~, t" 


1 

— I 5 — 


-1 


i 


-5 + 2sin^ 6 W 


u, c, t 




2 

-A" 




2 — ^ sm t/vK 


d, s, b 


2 


3 


-V 

2 


-| + f sin 2 6^ 



Table 1: Coefficients in Equation (fTOj) . 



P 



P 




Figure 31: Feynman diagram for pp — > H + X. 




f 



Figure 32: Feynman diagram for iJ + —>//'. 
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13 Production of H + 

From the diagrams in Figures and we obtain 
a(pp - H+X) = T ( H+ ~> 9*0 



dXg 

1A Xa 



V) • /|,(^,m 2 ) + /£,(*„, to 2 ) • /f(^,m 2 ) 



where 

7h = 



and 



r(H + -> //') = r(#- -> /'/) 

= ygg^l^g . A l/2 



m 2 m 2 , 



x {A 2 [to# - (to/ + to/') 2 ] + B 2 [m 2 H - (m f - m f >) 2 ] } 
with N c = 3 for quarks, iV c = 1 for leptons, 

A — mf tan (3 + to/ cot /3 

B = rrifi tan (3 — m,f cot /3 
/ = u,c,t,v e ,v^,v T , and /' = d, s, 6, e", r".[7j 

14 Production of h°W + X 

Let us now consider the channel pp — > /i°W /+ X. We obtain: 



dyd{p T f 



^2 / <&a[/g (x a , m 2 ) /? (x 6 , m 2 ) 



At 

where 

f p _ fP f ? — fP — fP f p — f? 

_ \fsm T e y + m 2 - m 2 ^ 

•Eamin /= ~ > 

s — y'srriTe v 
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m T = (m 2 w +Pt) 2 , 
_ x a ^/sm T e~ y + to|q - mw 

S X a Xi)Sj 

2 _ A (s,m^ ,m^) sin 2 6> 
Pt ~ 4l ' 

[m^o + m 2 ^ - s - cos 6A 1/2 (s, m 2 ^ , m 2 ^ )] , 
[rn^o + m 2 ^ — s + cos6'A 1/2 (s,m 2 l o,m 2 y)] , 



cos 6* = 1 — 



A(s,m 2 ,m 2 y ) 
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y is the rapidity of W + and px is the transverse momentum of W + . From 
the Feynman diagrams of Figure El we obtain for //' — > h°W + : 



\Vff>\ 2 G 2 F {\C H+ \ 2 sA 



da 1 

x [mj, tan 2 (3 + m) cot 2 p] + m% \C W \ 2 [%sm 2 w + A sin 2 9] 



—2Ch+^(Cw) [m 2 , tan/3(sA + 2m 2 v u (s — m 2 h o) 
+2m^/ {2m 2 h o — t)) — m 2 cot P(sA + 2m 2 y t (s — m 2 h o) 
+2m 4 w (2m 2 H0 -u))} 
1 

2" M! 



+-m^A sin 2 6* 



171 jC j TTlj/C j/ 



+2C H+ s 



mlorn^ — r A sin 2 9 



1 

4' 



rrv 



-S[m}C} +m },C},] 
cot (3Cf m 2 ,t3Ji(3Cfi 



u 



+2Ch+s [m 2 , tanjSCf 



iu — m 



cot (3C f i] 



-2&(CW) 



-A sin 9 m w + - 1 — sm hQ m w + Asm w + Am 



w 



fn 2 Cf m 2 ,Cfi 



t 



u 



- 2$t(C w )m 2 f C f [-2m 4 w + t(s- 2m 2 w )] 



-2^(C w )m 2 f ,C f/ [-2m A w + u (s - 2m 2 w )} } 
where A stands for A(s,m 2 i0 ,m 2 v ), 

cos (P — a) 



C 



w 



s 

sin (/3 — a) (s 



m 



H 



rn 



w 



im w T w ) 



[s — m 



w i 



C 



cos a 



sina 



i7) 



(89) 
(90) 



sin P cos P 

For pp — > h°W~X interchange u <-> t. 

For the Standard Model we obtain the differential cross section (|8*7|) with 
h° replaced by the Standard Model higgs, Ch+ = 0, Cf = Cf> = —1, and 
sm(p-a) = 1 in JHS|). 



15 Numerical examples 

Two sensitive channels for the search of the Standard Model higgs are pp — > 
h°ZX and pp — > h°W X. The cross section for pp — > h°ZX off resonance in 
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f ar ,( a\ — i n 

I dll I ju J — 1 U 


I dill fJ I — 


00 


0.20E+1 


0.13E-1 


0.35E-1 


bb 


0.31E+2 


0.31E+0 


0.12E-1 


cc 


0.73E-7 


0.73E-5 


0.18E-3 


ss 


0.20E+0 


0.20E-2 


0.81E-4 


dd 


0.10E-1 


0.10E-3 


0.41E-5 


uu 


0.18E-9 


0.18E-7 


0.46E-6 



Table 2: Production cross section [pb] for pp —>■ A from the indicated par- 
tons. m A o = 200GeV/c 2 , ^s = 1960GeV/c 2 . 



partons 


tan(/3) = 100 


tan(/3) = 10 


tan(^) = 2 


99 


0.42E+0 


0.22E-2 


0.19E-1 


bb 


0.82E+1 


0.82E-1 


0.33E-2 


cc 


0.19E-7 


0.19E-5 


0.49E-4 


ss 


0.57E-1 


0.57E-3 


0.23E-4 


dd 


0.44E-2 


0.44E-4 


0.18E-5 


uu 


0.89E-10 


0.89E-8 


0.22E-6 



Table 3: Production cross section [pb] for pp —>■ A from the indicated par- 
tons. m A o = 250GeV/c 2 , ^s = 1960GeV/c 2 . 

the Doublet model differs from the Standard Model by a factor sin 2 (/? — a) 
(see Equation (J7UJ)) and it will be hard to obtain both m h o and tan(/3). We 
are therefore interested in the production of h°Z on resonance. In particular 
pp — > A followed by A — > h°Z — > bbl~l + where I — fi,e. A peak should be 
observed in the h°Z invariant mass. From Equation (jfil)) we obtain the cross 
sections listed in Tables 121 and El 

Let us now consider the decays of A . As an example we take m h o = 
120GeV/c 2 , m H o = 250GeV/c 2 , m H = 200GeV/c 2 and m A o = 250GeV/c 2 . 
The corresponding branching fractions are listed in Table EJ From Tables 
H2 and E] we obtain a production cross section times branching fraction for 
the process pp -> A -> h°Z of 0.018pb for tan(/?) = 2, and 0.0045pb for 
tan(/?) = 10. 

From Equations ()71|) and (f?3*|) we obtain the production cross sections for 
pp — > H + X shown in Table El 

Other channels of experimental interest are the production of 3 or more 
6-jets as in Figure 1221 Some numerical calculations using the CompHEP 
program 8 are presented in Table |H1 
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partons 


tan(/3) = 100 


tan(/3) = 10 


tan(/3) = 2 


A^gg 


3.0E-4 


1.5E-4 


2.0E-3 


A^bb 


1.0E+0 


9.4E-1 


5.9E-2 


A—>cc 


8.0E-10 


7.5E-6 


2.9E-4 


A^ ss 


8.0E-4 


7.5E-4 


4.7E-5 


A — > Zh° 


6.1E-6 


5.5E-2 


9.4E-1 


A—>Z~j 


1.9E-8 


6.0E-9 


1.2E-6 


A — > 77 


1.2E-7 


8.1E-8 


7.5E-6 



Table 4: Branching fractions for A assuming uih = 200GeV/c 2 , m#o = 
250GeV/c 2 , m A o = 250GeV/c 2 and m^o = 120GeV/c 2 . 



partons 


tan(/3) = 100 


tan(/3) = 10 


tan(/3) = 2 


ud 


0.82E-2 


0.82E-4 


0.34E-5 


us 


0.66E-1 


0.66E-3 


0.26E-4 


ub 


0.89E-2 


0.89E-4 


0.36E-5 


cs 


0.31E-1 


0.32E-3 


0.91E-4 


cb 


0.24E-1 


0.24E-3 


0.96E-5 



Table 5: Production cross section [pb] for pp — > iJ + X from the indicated 
partons. m H o = 250GeV/c 2 , y/s = 1960GeV/c 2 . 



process 


tan(/3) = 100 


tan(/3) = 50 


tan(/3) = 2 


b + g ^ b + h° 


0.021 


0.021 


0.011 


u + u^b + b + h 


0.002 


0.001 


0.0004 


d + d^b + b + h° 


0.0005 


0.0005 


0.0001 


g + g -> 6 + 6 + /i° 


0.015 


0.015 


0.008 



Table 6: Production cross section [pb] for pp — > &/i°X from the indicated 
processes. m h o = 120GeV/c 2 , m A o = 250GeV/c 2 , ^1 = 1960GeV/c 2 . 
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Figure 33: Some Feynman diagrams for the production of three or more 
6-jets. 
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16 Running coupling constants and Grand Uni- 
fication 

The coupling constants of the Two Higgs Doublet Model of type II are g s {fJ>) 
for SU(3), g(n) for SU(2), and g'(fj) for U(l). These coupling constants 
depend on the energy scale \i as follows: 



+ 



(91) 



+ 



g 2 (ji) g 2 (m x ) 8n 



22 4 1 
-— + -n F + -n s 
3 3 6 



, ( m-x 
In — 

V u 



(92) 



1 



20 1 
— n F + -n s 
9 o 



In — 



g*Qi) g'\m x ) ' W \ 9 "* ' 6'"V ~ W ^ 

where rip is the number of families of quarks and leptons, and ns is the 
number of higgs doublets. For the Two Higgs Doublet Model of type II con- 
sidered in this article, rip = 3 and ns = 2. In terms of the elementary elec- 
tric charge and the Weinberg angle, g(mz) = e{m z )/ sva.9wi.rnz), g'ijnz) = 
e(m z )/ cos6 w {rn z ). The fine structure constant is a(mz) = e 2 (m^)/(47r). 

Let us now assume that a Grand Unified Theory (GUT) breaks its sym- 
metry to SU(3)xSU(2)xU(l) at the energy scale m x . At this scale we take 



and obtain 



g 2 s {m x ) 



g 2 {m x ) 



3 9 



'2, 



m x ), 



sin 2 6 



11 + iris + g (22 - in 8 ) 



w 



In 



m x 
m z 
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66 + n s 
24tt 2 1 " 



in 



e 2 66 + ns ' 



(94) 



(95) 



(96) 



with all running couplings evaluated at mz- 

The corresponding equations of the Minimum Supersymmetry ModeljOj 

are 

18 + 3n s + $ (60 - 2n s ) , 
sin 2 W = T7 ^ 7 , (97) 



108 + 6n s 



In ( — 



57T 



1 8e 2 
3£f_ 

18 + n 5 



(9* 
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Doublet Model 


MSSM 




sin 2 6 w {m z ) 


m x 


sin 2 6 w {m z ) 


m x 





0.2037 


1.0 ■ 10 15 


0.2037 


8.0 • 10 17 


2 


0.2118 


4.2 ■ 10 14 


0.2311 


2.0 • 10 16 


4 


0.2194 


1.8 ■ 10 14 


0.2536 


1.0 ■ 10 15 


6 


0.2266 


8.3 • 10 13 


0.2722 


8.3 • 10 13 


8 


0.2334 


3.9 • 10 13 


0.2880 


1.0 • 10 13 



Table 7: Predicted sin 2 9w{ m z) an d m x for the Two Higgs Doublet Model 
of type II, and the Minimum Supersymmetric Model as a function of the 
number of doublets ns- 

Some numerical results are presented in Table [H From the Table we 
conclude that the Two Higgs Doublet Model of type II is in disagreement 
with the measured value of sin 2 ^^y(m^), and with the non-observation of 
proton decay (m x is too low). Raising the number of doublets to ~ 7 would 
bring sin 2 #iy( m z) into agreement with observations, but m x is still too low. 
The MSSM with ns = 2 (which includes the Two Higgs Doublet Model of 
type II) is in agreement with both the observed sin 2 #^(mz), and with the 
non-observation of proton decay. 

17 Conclusions 

One of the major efforts at the Fermilab Tevatron in Run II, and at the future 
LHC, is the search for the Standard Model higgs Hsm- The four channels with 
largest production cross section arejO] gg — > hsM, qq' — > hsntW, qq — ► HsmZ 
and qq — > hsmqq. The decay modes of hsM with largest branching fractionjOj 
are bb for m h < 137GeV and W+W~ for m h > 137GeV. 

The search for the Standard Model higgs will also constrain or discover 
particles of the Two Higgs Doublet Model of type II. 

The most interesting production channels are gg -> h°, H°, A on mass 
shell, and qq, gg h°Z and qq' — > h°W ± in the continuum (tho there may 
be peaks at m^o). The most interesting decays are h°, H°, A bb-jets and 
t + t~, and, if above threshold, H° — > ZZ, W + W~ and h°h°. The following 
final states should be compared with the Standard Model cross section: bbZ, 
bbW ± , t + t-Z, t + t~W ± , bb, t + t-, ZZ, W + W~, 3 and 4 6-jets, 2r+ + 2r~, 
bbr+r-, ZW + W~, 3Z, ZZW ± and 31^ Mass peaks should be searched in 
the following channels: Zbb, ZZ, ZZZ, bb, 46-jets and, just in case, Zj. 

We have discussed the masses of the higgs particles in the Two Higgs 
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Doublet Model of type II, and have calculated several relevant production 
and decay rates. We have also discussed running coupling constants and 
Grand Unification. If the Two Higgs Doublet Model of type II is part of a 
Grand Unified Theory, then it does not agree with the observed sin 2 9w nor 
with the non-observation of proton decay. The MSSM with ns = 2 (which 
includes the Two Higgs Doublet Model of type II) is in agreement with both 
the observed sin 2 9w(fnz), and with the non-observation of proton decay. 
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